
use Data/Andrea/sample_income_discountoc_angle.dta, clear

 keep if read_date > = td(01jan2013) // same picture if restrict to > jan2013, feb2013, mar2013


collapse (mean) mean_ghi*  ghi Melbourne (firstnm) FIT_t FIT_p streetangle  , by(account_number)

 
gen year = 2007 
replace year = 2009 if FIT_t=="P"
replace year = 2011 if FIT_t=="T"
replace year = 2013 if FIT_t=="F"


label define year_lab 2007 "1:1 FIT" 2009 "60 FIT" 2011 "25 FIT" 2013 "8 FIT"
label values year year_lab







graph box ghi, over(year) ytitle("Mean hourly irradiance") box(1, color("178 205 248"))  ///
 box(3, fcolor("255 250 183")) box(4, fcolor(gs9)) nooutsides scheme(s2mono) graphregion(color(white))
 
 graph export Results/Andrea/mean_irradiance.eps, replace
 graph export Results/Andrea/mean_irradiance.png, replace

 
gen not_NS= streetangle~=3
label variable not_NS "Street angle"



 graph drop _all
 graph box ghi, over(year) ytitle("Mean hourly irradiance", size(large))  ///
  scheme(s2mono) graphregion(color(white)) name(ghi) title("Panel A. Sunshine availability", position(11) size(large) ring(7))
 graph export Results/Andrea/mean_irradiance2013.pdf, replace

 
 
 
 
graph bar not_NS, over(year) ytitle("Street angle", size(large))  ///
  nooutsides scheme(s2mono) graphregion(color(white)) name(angle) title("Panel B. Street angle", position(11) size(large) ring(7))
 graph export Results/Andrea/streetangle.pdf, replace

 
 
 graph combine ghi angle 
 
 graph export Results/Andrea/ghi_angle.png, replace
 
 
